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Abstract: 

We study power corrections to the differential thrust, heavy jet mass and C— parameter 
distributions in the two-jet kinematical region. We argue that away from the end-point 
region, e 3> Aqcd/<5, the leading l/Q— power corrections are parameterized by a single 
nonperturbative scale while for e ~ Aqcd/Q one encounters a novel regime in which 
power corrections of the form l/(Qe) n have to be taken into account for arbitrary n. 
These nonperturbative corrections can be resummed and factor out into a universal 
nonperturbative distribution, the shape function, and the differential event shape distri- 
butions are given by convolution of the shape function with perturbative cross-sections. 
Choosing a simple ansatz for the shape function we demonstrate a good agreement of 
the obtained QCD predictions for the distributions and their lowest moments with the 
existing data over a wide energy interval. 
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1. Introduction 



Analysis of hadronization effects to the final states in e + e~— annihilation has became the 
subject of active QCD studies jj]]. There exist infrared and collinear safe event shape 
variables for which perturbative QCD can be applied at large center-of-mass energies 
s = Q 2 to calculate their differential distributions and mean values as series in a s (Q). It 
has been observed many years ago || that for some shape variables like thrust, t, and 
heavy jet mass, p, perturbative QCD predictions deviate from the data by corrections 
suppressed by powers of the large energy scale l/Q p , with the exponent p depending 
on the variable and p = 1 for t— and p— variables. Such hadronization corrections were 
measured experimentally JTJ] over a wide energy interval 14 < y/s/GeV < 189 and were 
found to have a different form for the differential event shape distributions, da/de, as 
compared to their mean values, (e) = a^l J deeda/de. For the mean value (e) the 
leading power correction is parameterized by a nonperturbative scale X p of dimension p, 
while hadronization corrections to the differential distribution are described by a function 
/hadr(Q, e) depending on both the shape variable and the center-of-mass energy 

(e) = (e) PT + \ P /Q P , = ^ + / hadr (Q, e) (1.1) 

with e denoting a general event shape variable (e = t,p,C ,...) and the subscript PT 
referring to perturbative contribution, (e) PT = J deeda PT /de. Obviously, the hadroniza- 
tion corrections to the differential distributions have a richer structure then those to the 
mean values. For instance, nonperturbative scales X p parameterizing power corrections 
to (e) are defined by the moment J de e/hadr(Q, e). 

Power corrections in (|1 . 1| ) are associated with hadronization effects in e + e~— final 
states and, as a consequence, the magnitude of the scales A p and the function fhadr{Q, e) 
cannot be calculated within perturbative QCD approach. However it was recognized 
some time ago || ||, [| [J, that analysis of infrared renormalon ambiguities of perturba- 
tive QCD series suggests the value of dimensionless exponents p as well as the dependence 
of the function /hadr(<3, e) on the large scale Q. Namely, perturbative QCD series gen- 



erate power corrections of the form (1.1) through IR renormalons contribution but fail 
to predict uniquely their values - it is only the sum of perturbative and nonperturbative 
contributions that becomes well-defined 0. To give a meaning to the perturbative se- 
ries in ( |1.1| ) one has to regularize IR renormalon singularities. This can be done in two 
different ways: one can specify a particular prescription for integrating IR renormalon 
singularities like principal value prescription ||. Alternatively, one can avoid IR renor- 
malon ambiguities by introducing an explicit IR cut-off p on momenta of soft particles 
in perturbative expressions. In this case, one can either impose a "hard" IR cut-off on 
momenta of soft particles in the Feynman integrals, k± > p, || or replace QCD coupling 
constant by a effective IR finite coupling constant which coincides with a s (kj_) at large 
scale k± and deviates from it at k± < p [Q, [5|. Following each of these ways, one speci- 
fies perturbative (//—dependent) contribution to fll.lj ) including perturbatively induced 
power corrections. Still, there exists a genuine nonperturbative contribution to the event 
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shapes coming from the QCD dynamics at scales below fi. This contribution cannot be 
determined from the analysis of perturbative QCD series while its magnitude depends 
on the choice of the IR regularization and, as a consequence, on the IR cut-off \i. 

For some hadronic observables like mean values of the event shapes and their dif- 
ferential distributions away from the end-point region, leading nonperturbative power 
corrections can be parameterized using different IR renormalon inspired phenomenolog- 
ical models 0, ^ f|, [7J. Their predictions agree well with the experimental data and the 
extracted values of phenomenological nonperturbative parameters exhibit approximate 
universality. Despite a phenomenological success of these models, it remains still unclear 
what is the physical meaning of new nonperturbative QCD scales and what is the origin 
of the universality property within QCD. In the present paper we address these prob- 
lems using the factorization properties of the event shape distributions established in [|Ttj . 
We shall argue that nonperturbative power corrections to the thrust, heavy jet mass and 
C— parameter distributions are described by the universal shape function which is a new 
nonperturbative QCD distribution measuring the energy flow in the two-jet final states 
in e + e~— annihilation. 

The paper is organized as follows. In Sect. 2 we discuss the general properties of 
power corrections to the event shape distributions in the end-point region. In Sect. 3 
we formulate the factorization procedure and define the shape function. In Sect. 4 we 
show that the differential event shape distributions are given by the convolution of the 
resummed perturbative cross-sections with universal shape function. Choosing a simple 
ansatz for this function we compare QCD predictions with the existing data. In Sect. 5 
we apply the obtained expressions to calculate the power corrections to the first two 
moments of the distributions. Concluding remarks are given in Sect. 6. 



2. Event shape distributions 

In this paper we shall consider three event shape variables: thrust T, heavy jet mass p 



and C— parameter. They are defined in the standard way as [ITT 



rp Ek \Pk-fh\ (Ml Ml\ 

T = max — — , p = max — f , — , (2.1) 



J2 k \Pk\ ' \Q 2 ' Q 



where M\ and M\ denote the total invariant masses flowing into the right and left hemi- 
spheres with respect to the plane orthogonal to the thrust axis ut- The C— parameter 
is given by 

C = 3(0^ + 2 3 + 3 0!) (2.2) 
with 6j being eigenvalues of space-like part of the energy-momentum tensor Q alS = 



3 I 



T, k PkPkl\Pk\IT, j \p. 

Introducing the new variable t = 1 — T one notices that thus defined event shapes e = 
(t , p ,C) have a number of common features. Lowest order perturbative QCD calculation 
leads in all three cases to the following expression for the differential distribution for e > 
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de 



<*s{Q) 

2tt 



A e (e)9(e max - e) + 



<*s(Q) 
2tt 



B e (e) + 0(a 3 s ) 



(2.3) 



where v4 e and _B e are known coefficient functions and normalization is chosen as J de 
1. Lowest order correction A e gets contribution only from the three-particle final state 
which populates the kinematic region < e < e max with t max = p max = 1/3 and 
Cmax = 3/4. Away from the end-point region, e A QCD /Q, the perturbative expansion 
Q2.3| ) is well-defined and it describes the final states consisting of particles with relative 
transverse momentum that scales at large center-of-mass energy as ~ Q. As e approaches 
the three-particle upper limit, e = e max , A t and A„ vanish while Ac takes a finite value 
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A t (l/3) = A p (l/3) 
256 



A C (C) 



-71 



243 



C--)+0((C-3/4) 2 ) 



(2.4) 



For e = (t ,p,C) — > A QCD /<5 the final states consist of two narrow jets with invariant 



mass M| L ~ 



A 



QCD 



as 11 



mm 



Q. Examining 

Tl 



273|) one finds that A e diverges in the end-point region 



AJe) 



4C f 



+ C(lne) 



(2.5) 



with to = po = 1 and C = 6. Similar Sudakov-type corrections appear to higher orders, 



a 



N \n 2N ~ n e/e with n > 0, and need to be resummed |TB|, [T3j]. They originate from the 
effects of collinear splitting of quarks and gluons inside two narrow energetic jets and 
their interaction with surrounding cloud of soft gluons. The underlying QCD dynamics 
depends on two infrared scales, Qe and Q 2 e, such that 1/Q l/(Qe l l 2 ) l/(Qe). 
The smallest scale Qe sets up the typical energy carried by soft gluons, while the scale 
Q\fe defines the transverse momenta of the jets, k\ = Q 2 e. Applying the standard IR 
renormalon analysis and examining sensitivity of perturbative expressions with respect 
to emission of particles on each of these scales, that is soft gluons with energy ~ Qe and 
collinear particles with the transverse momentum ~ Q 2 e, one finds that nonperturbative 
corrections to the differential distribution appear suppressed by powers of both scales. 
Then, in the end-point region, e = C(A QCD /Q), we may use the fact that Qe = 0(A QCD ) 
and expand the differential distribution in powers of larger scale Q 2 e. Keeping only the 
leading term of the expansion one gets 



1 da 

a tot de 



cr (a s (Q),lne, 



Qe 



O 



1 



Q 2 e 



[2.6) 



where do resums perturbative corrections in a s (Q) as well as power corrections on the 
smallest scale Qe 



0o 



G?<T pl 

de 



E 



A* 



{ {Qe)' 



: S fc (a s (Q),lne) 



(2.7) 
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Here, are dimensionless perturbative coefficient functions and A& are some nonper- 
turbative scales, depending, in general, on the choice of the event shape variable. Using 
(p.7|) we notice that the power corrections have a different form for e 3> A QCD /Q and 

e ~ A qcd/Q- 

For e away from the end-point region, e ^> A QCD /Q, one may keep in (|2.7f) only the 
first term 

= io^ + v / i \ 

atot de de Qe \{Q e ) J 

The coefficient function Si can be found using the well-known property |4], [15|] that the 
leading 1/Q— power correction to the differential distribution ( |2.8| ) is generated by a shift 
of perturbative spectrum, e — ► e — X±/Q. This leads to 



Si(a s (Q),lne) 



d 



dine 



dov 



de 



(2.9) 



Then, it follows from (|2.8| ) that for e 3> A QCD /Q the leading power corrections to the 
differential distributions have a rather simple structure: they are parameterized by a 
single nonperturbative scale Ai. The same scale determines 1/Q— power correction to 
the mean value (e). The QCD predictions Q2.8| ) are in a good agreement with the 
experimental data and the value of Ai has been fitted for different shape variables [Ij]. 
It is worthwhile to note that in the performed analysis of power corrections to the 
differential distributions |3J the fitting range of event shape variables was restricted to 
the region e ^> A QCD /Q. At the same time, as we shall argue below, it is in the region 
e ~ A QCD /<5 where a novel QCD regime is realized and the structure of hadronization 
corrections is drastically changed. 

For e ~ A QCD /Q one finds that all terms in ( |2.7| ) become equally important and, 
therefore, need to be resummed to all orders in l/(Qe). The resummation is based on 
the remarkable factorization properties of the differential distributions. As was shown 
in [|K| , the nonperturbative corrections to the leading asymptotic term ctq are factorized 



out into nonperturbative distribution function, the so-called shape function. The general 
factorized expression for differential distribution looks like |TU| 



1 da 

CTtot de 



def e (e 



dcr pT (e - 
de 



+ o 



Q 2 e 



(2.10) 



Its explicit form depends on the choice of the event shape variable e = (t,p,C) and will 
be given in Sect. 4 (see Eqs. ( |4.1| ), ( |4.2j ) and (|4.3|) ). For e ^> Aqqd/Q one can expand 
the r.h.s. of ( 2.10Q in powers of 1/Q to reproduce the expansion (|2.7|) with 



A, 



dee n f(e) 



S„(a s (Q),lne) 



-e) n d n 



ni 



de n 



da^ 



de 



(2.11) 



Expression fl2.10|) has a simple physical interpretation - nonperturbative corrections 
increase invariant masses of jets and effectively shift perturbative spectrum towards 
larger values of the shape variables with the weight given by nonperturbative distribution 

f e (e). 
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3. Factorization and Shape functions 



Factorization relations (|2.10 ) take a simple form for the radiator functions R(e) defined 
as [|13[ 

1 da 



R{e)= / de >_ = {B{e _ e (jv))) . (3.1) 
' cr tot de' 



o 



Here, (...) denotes averaging over all possible final states in e + e~— annihilation with the 
weight given by the differential distribution 1/ 'a t otd<J / 'de and e(N) denotes the value of 
the event shape variable e for a given final state \N). Calculating R(e) in perturbation 
theory one finds 

R PT (e) = 1 - *M F~ de'A e (e') + 0(a 2 s (Q)) . (3.2) 
2tt J e 

Close to the two-jet region, e — ► 0, perturbative expressions for i?(e) involve Su- 
dakov logs af vf ln 2Ar-n e with n > 0. In the case of event shapes under consideration, 
e = (t , p ,C), these corrections can be systematically resummed to next-to-leading log- 
arithmic (NLL) order and matched into exact two-loop perturbative expressions ( |3.2j ) 
fl3| , |i4f . To this accuracy the weights e(N) can be expressed in terms of the total 
invariant masses M| and M| of two jets flowing into the right and left hemispheres, re- 
spectively. Moreover, the t— and C— parameters depend only on the sum of two masses 
and the corresponding perturbative radiation functions can be expressed to the NLL 
approximation as [TJJ] 

/ / m 2 + /If 2 \ \ 
i?f T (e) = (B(e - t(N))) PT = (o[e- R ^ 2 L j j (3.3) 

R-J{e) = (6(e - C(N))) PT = (<>(*- 6 ^^ ) ) = ^ PT (e/6) , 

where the subscript PT indicates that the final states in e + e~— annihilation are generat- 
ing by perturbative branching of outgoing quark and antiquark. The radiator function 
for the p— parameter depends separately on the masses of two jets. Taking into ac- 
count that perturbative evolution of two jets is independent on each other to the NLL 
approximation one gets |T3| 

U 7(e) = W e - ^))> PI = /,(.- f ) ) (, L - f ) ) (3.4) 



The perturbative expressions ( |3.3| ) and ( |3.4| ) are valid in the two-jet kinematical region 



A QCD /Q -C e < e max except the end-point region e ~ A QCD /Q, in which the energy of 
emitted soft particles scales as k± ~ eQ ~ A QCD and perturbation theory is expected to 
fail. 

Calculating the radiator functions R e {e) one has to combine together perturbative 
and nonperturbative corrections. In the case of inclusive distributions, like deep inelas- 
tic structure functions and Drell-Yan distributions, this can be achieved by applying the 
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factorisation theorems. They allow to separate short- distance dynamics into perturba- 
tively calculable coefficient functions and absorb large-distance corrections into universal 
nonperturbative distributions. Specific feature of the differential event-shape distribu- 
tions is that they are not inclusive quantities but rather weighted cross-sections and, as 
a consequence, the standard methods are not applicable in this case. 



It turns out [JXO | that IR factorization still holds for the leading term o" in the expan- 
sion of the event-distributions ( |2.6| ) in the end-point region e ~ Aq CD /Q. Its origin has a 
simple physical interpretation. In end-point region, the final state in e + e~— annihilation 
consists of two narrow jets surrounding by a cloud of soft gluons. Nonperturbative cor- 
rections ~ 1/ (Q 2 e) and ~ l/(Qe) are associated with emission of collinear particles with 
the transverse momenta k\ ~ Q 2 t and soft particles on the energy scale k± ~ Qe, re- 
spectively. Neglecting power corrections to ( |2.6| ) on a larger scale, ~ l/(Q 2 e), we may 
restrict analysis to soft particles only. Since soft particle cannot resolve the internal 
structure of narrow jets of transverse size k\ ~ Q 2 e, we may effectively replace two jets 
by a pair of energetic quark and antiquark moving back-to-back with the energy ~ Q/2. 
The internal dynamics of two jets is governed by perturbative branching of quark and 
antiquark while effects of their interaction with soft gluons can be factorized out into the 
eikonal phase W+iyl with W + and W- being the eikonal phases of quark and antiquark, 
respectively. They are given by Wilson lines W± = Pexp(i J °° dsn±A(sn±)) in which 
soft gluon field A^x) is integrated along the light-like directions n± defined by the mo- 
menta of two outgoing jets. In the end-point region, collinear and soft particles provide 
additive contributions to the shape variables e = (t , p ,C). As a consequence, the radia- 
tor functions are given in all three cases by a convolution of perturbative radiators R PT 
and the same universal nonperturbative distribution /(er, El) describing the energy flow 
into the right and left hemispheres in the final state, Er and El, respectively, created 
by nonperturbative soft gluon radiation. The nonperturbative distribution /(e^el) is 
defined as follows 111 



f(s R ,s L ) = J2mW + Wl\N)\ 2 6(E R - (k R n + ))5(E L - (k L nJ)) . (3.5) 

N 

Here, sum goes over all possible soft gluon final states \N) with k R and ki being the 
total momentum of soft particles moving into right and left hemispheres, respectively. 
The quantities {k R n + ) and {k L nJ) define the projection of the soft gluon momenta onto 
the directions of two jets, n^. = (1, 0j_, ±1), propagating into the same hemisphere. 

Finally, the factorized expressions for the radiator function for the t— and C— variables 
look like 

(3.6) 



Rt(e) = £ dE f t (e) Rr (e - £) 



'£ eQ ( 3tt £ 

UrU)= I dEf t {E)R-j[e---) (3.7) 
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with nonperturbative distribution f t (e) defined as 



ft{e) = J de R J de L f(e R ,e L ) 5(e - e R - e L ) = de' f(e - e',e') . (3.8) 
In the case of the p— variable, 

R P {e) = jT de R J* Q de L f(e R , e L )E? (e - ^ R p / (e - E ± 



(3.9) 

with R p p T (e) = [R p j T (e)} 2 and R p /{e) = (9 (e - M R /Q 2 )) pT being a single jet radiator 
function. We would like to stress that Eqs. ( |3.6|) , ( |3.7|) and (|3.9|) hold in the region 



Aq CD /Q 2 < e < e max . They resum all power corrections of the form l/(Qe) n and 
are valid up to corrections ~ l/(Q 2 e). According to (|3.6|) , (|3.7| ) and ([3.9|) , the power 
corrections have a different form for p and e = (t, C) variables. In the latter case, 
the radiator function depends on an overall energy flowing into both hemispheres and 
described by the integrated distribution (|3.8|). 



Nonperturbative corrections to the radiator functions ( |3.6|) , (|3.7|) and (|3.9| ) are gov- 
erned by the universal shape function f(s R , e£). This function is different from the 
well-known inclusive QCD distributions and its operator definition was given in |JL(]|| . 
Using ( |3.5| ) it is straightforward to show that f(s R , El) is a symmetric function of its 
arguments, it does not depend on the center-of-mass energy Q and is normalized as 



d 



dQ' 



f(£R,£L)=0, J de R J de L f(e R , e L ) = 1 , (3.10) 



where the last relation follows from unitarity of the eikonal phase W + W}_. The matrix 
element entering (|3.5|) does not depend on any kinematical scale and, as a consequence, 



the momenta of soft gluons contributing to (|3.5|) are not restricted from above. To sep- 



arate the region of small gluon momenta one has to introduce the factorisation scale 
p. Then, the shape function describes the contribution of gluons with k± < p, while 
the contribution of gluons with k± > p, is absorbed into perturbative radiator func- 
tion R(e). In this way, both nonperturbative shape function and perturbative radiator 
become //—dependent while this dependence cancel in their convolution ( |3.6|) , ( |3.7|) and 
(^.9[). Since the p— dependence of radiator function R PT can be calculated perturbatively, 
the above condition allows to obtain the evolution equations on the nonperturbative dis- 
tributions fllQ"| . Clearly, there exists an ambiguity in implementing IR cut-off inside 
perturbative expressions. Different prescriptions correspond to different ways of regular- 
izing IR renormalon singularities and therefore lead to the different expressions for the 



nonperturbative distributions. In what follows we shall impose a "hard" IR cut-off |T0 



on gluon momenta inside the perturbative radiator functions entering (|3.6|), ( p.7| ) and 
as [0 

R PT (e) -> R PT {e-p) = o(e- ^)iC» + -«)CW«) ■ (3.H) 
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Throughout the paper we shall substitute i?^ L (e) by its perturbative expression re- 
summed to the NLL accuracy and matched into two-loop explicit expressions within the 
modified In R— matching scheme [ 13| . Thus defined radiator function (|3.11|) depends on 



two scales, Aqcd and IR cut-off ll, that we choose as 

Aqcd = /U = 0.25 GeV. (3.12) 



Within the prescription ( ft.llj ), the "regularized" perturbative spectrum da PT (e; jj)de 



dR PT (e; ll) / de coincides with the In R— matched perturbative distribution dR (e) / de 
for li/Q < e < e max and it vanishes inside the nonperturbative "window" < e < 
li/Q. Choosing the value of \i in one has to be sure that the end-point of the 

perturbative distribution, e = li/Q, belongs to applicability range of the NLL resummed 
radiator function R pT (e) 2j3 a s (Q 2 ) lne < 1. Despite the fact that the perturbative 
spectrum is well defined at e = li/Q we do not expect that it provides a reasonable 
description of the physical distribution in the end-point region. Indeed, it is in this 
region that nonperturbative power corrections become dominant. 



4. Differential distributions 



Differentiating the radiator functions (|3.6| ) and (|3.7|) we obtain the following expressions 
for the differential t— distribution 



1 da t 

0tot de 

and C— distribution 



Q/(Qe;//)i?r(0;//)+ / def t {e;n) 

Jo 



da f T (e- e/Q;ll) 
de 



(4.1) 



1 da, 



c 



(Jtot de 



R p J(0;ll) + 



zQe 



deftfafi)- 



da p J 



2 Q> V 



de 



(4.2) 



Here, we indicated explicitly the dependence of nonperturbative shape function and 
perturbative distributions on the factorization scale ll. Two terms entering the r.h.s. 
of ( |4.1|) and ( |4.2| ) have the following interpretation. Since the shape function ft{s) 
rapidly vanishes for large e, the first term contributes inside the nonperturbative window 
< e < li/Q. In this region the emission of perturbative real soft gluons is suppressed 
due to cut-off imposed on soft gluon momenta k^> ll and the shape of the distribution is 
governed entirely by nonperturbative function ft(e). Additional Sudakov factor i? PT (0; ll) 
takes into account the contribution of virtual soft gluons with ll < k± < Q and it rapidly 
vanishes as ll decreases. The second term in ( |4.1| ) and ( |4.2j ) defines the spectrum inside 
the perturbative window li/Q < e < e max . In this region, nonperturbative corrections 



smear the perturbative spectrum over the interval Ae 
For the heavy mass distribution one finds 



^qcd/ Q- 



1 do p 

de 



Qf p (eQ,eQ; fJ .)RY(0; f i)+ / de fJe, eQ; ll) 



eQ 



da p j T (e - s/Q; ll) 
de 



(4.3) 
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where da^/de is single jet distribution resummed to the NLL order and defined by 
the radiator function ( |3.9| ), daj T /de = dR v j T (e)/de. The heavy mass nonperturbative 
distribution is given by 



/ p (£,eQ)= 2^ dE'f(e,e')Ryi E ' 



Comparing (|4.3| ) with (4J.) and (^2) we notice that the factorized expressions for 



the differential t—, C— and p— distributions have a similar form but the structure of 



power corrections is different in the case of the heavy mass. In distinction with ( |3.8| ), 
the heavy mass nonperturbative function f p depends on the shape variable, e, and 
the center-of-mass energy, Q. This dependence is controlled by perturbative radia- 
tor function and has the following interpretation. In the two-jet limit, the invariant 
mass of each jet is given by the sum of perturbative and nonperturbative contribu- 
tions, M\ = M^ PT + erQ and M| = PT + elQ- Perturbative radiation leads 
to Mft PT /Q 2 ~ Ml PT /Q 2 ~ 0(a s (Q)), while nonperturbative contribution scales as 
£r ~ £l ~ C(Aqcd)- In contrast with the t— variable, which depends on the sum of 
both masses and therefore is additive with respect to perturbative and nonperturbative 
contributions, the p— parameter is defined by the largest mass for which the "additiv- 
ity" property is lost. Namely, comparing invariant masses flowing into two hemispheres 
one encounters a situation when masses of two perturbative jets are of the same order, 
M r/l, pt = ®(Q 2 )i while their difference is much smaller \M 2 LPT - M 2 RPT \ = C(QA QCD ).0 
In this case, nonperturbative correction to the difference of the jet masses becomes 
comparable with the perturbative contribution |M| PT — M RPT \ ~ Q\sl — £r\, and 
therefore it can invert the perturbative hierarchy of jet masses, M RPT < M| PT , into 
M\ > Ml, for instance. Expression ( |4.3|) takes into account this effect through the 



induced Q— dependence of the nonperturbative function (4.4). 

According to (flip, (fO|) and ( |4.3| ) the nonperturbative corrections to the t—, C— 
and p— distributions involve two different nonperturbative functions. They are related 
however to the same universal nonperturbative shape function ( |3.5|) describing the energy 
flow into two hemispheres in the final state. We recall, that /(^h^l! ^) depends on the 
cut-off p imposed on the maximal momenta of soft particles but it is independent on the 
center-of-mass energy Q. By the definition, f [er. El] p) distinguishes between particles 
propagating into right and left hemispheres in the final state and therefore it is not 
completely inclusive with respect to partonic final states. Namely, it takes into account 
that quarks and gluons produced at short distances ~ 1/Q and moving into one of 
the hemispheres will eventually decay at large distances ~ 1/Aqcd and their remnants 
could flow into opposite hemispheres.^ This implies that, firstly, in contrast with the 
well-known inclusive QCD distributions, the shape function f(£R,£L]p) is not related 
to the short distance QCD dynamics and, in particular, its moments can not be related 
to hadronic matrix elements of local composite operators. Indeed, according to the 



2 To see that this configuration is not rare it is enough to notice that it corresponds to the vicinity 
of peak of the perturbative distribution over the difference of the jet masses \M\ — M^\/Q 2 H] , 
3 Similar effect has been studied using the IR renormalon approach in Jig. 
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operator definition proposed in [pi, the shape functions are defined in terms of the so- 
called "maximally nonlocal" QCD operators [|17], |18| . Secondly, non-inclusive corrections 
to the shape function describe a "cross-talk" between two hemispheres in the final state 
leading to correlations between er and sl- As a consequence, the shape function is not 
factorizable into the product of functions depending on the energy flowing into separate 
hemispheres 

f(e R , e L ) = j i nc i(e R )j \ nc i{e L ) + Sf non ^ mc i(e R , e L ) . (4.5) 

One should notice that similar property holds for perturbative Sudakov resummed radia- 
tor function (3^). However, one finds that there the factorization holds to the NLL accu- 
racy, Eq. (|3.4j ), and non- inclusive corrections first appear at the NNLL level ~ a 2 (a s L) N . 

In what follows we shall rely on a particular ansatz for the shape function f(e R ,ei) 
which agrees with general properties of nonperturbative QCD distributions and has been 
used in previous studies of power corrections to the thrust distributions |IIJ. Namely, 
one expects that for small values of e R and El the shape function should vanish as a 
power of the energy. Similarly, f(e R ,EL) should rapidly vanish as Er or el becomes 
large. Taking into account these properties together with ( |4.5|) one chooses the following 
expression 



Af(a,b) (ERELy- 1 ( e 2 r + e 2 L + 2be R e L \ 
^ = -^T- {-jr) exp { T2 ) 



(4.6) 



It depends on two dimensionless parameters a and b and the scale A. The factor Af(a, b) 
is fixed by normalization condition (|3.10|) . 

The free parameters, a, b and A, have the following meaning. The exponent a deter- 
mines how fast the shape function vanishes at the origin. The scale A sets up the typical 
energy of soft radiation. The parameter b controls the non-inclusive contribution to the 
shape function and its possible values are restricted as b > — 1 in order for the shape func- 
tion fl4.6| ) to be normalizable. Non-inclusive corrections vanish at b = 0, 5/ non -mci = 
in Q4.5|) . For b ^> 1, the shape function enhances the regions of the phase space £r^> el 
and Er -C £l, in which most of the energy flows into one of the hemispheres. For 6^—1 
the energies are of the same order, Er ~ El, and strongly correlated to each other. We 
expect that non-inclusive corrections to the shape function should be important and the 
configurations in which energy flows mostly into one of the hemispheres to be suppressed. 
This suggests that the possible values of the b— parameter should lie within the interval 
-1 < b < 0. 

The parameters a, b and A depend on the factorization scale fi and are independent 
on the center-of-mass energy Q as well as the choice of the shape variable e = (t,p,C). 
This allows to fit their values by comparing the event shape distributions, Eqs. ([4.1| ), 
( |4.2| ) and ([4.3|) , with the most precise available experimental data at Q = Mz- Following 
this procedure we found that the fit to the heavy jet mass distribution is more sensitive 
to the choice of the parameters (especially to the non-inclusiveness parameter b) then the 
one to the thrust and the C— parameter. Then, fitting the heavy jet mass distribution 
at Q = Mz as shown in Fig. |I]a we obtain 

a = 2, b = -0.4, A = 0.55GeV. (4.7) 
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Using these values we compare the QCD predictions for the C— parameter distribution 
at Q = Mz with and without nonperturbative corrections included as shown in Fig. |l|b. 
Similar plot for the thrust distribution can be found in ||10|| . We observe that the dif- 
ferential distributions ( |4.3|) and ( f4.2|) combined with the shape function, Eqs. (|4.6| ) and 
( flT]) , correctly describe the data throughout the interval < e < e max including the 
end-point region e = C(Aq C d/Q)- In addition, the p— parameter distribution turns out 
to be very sensitive to the choice of the b— parameter. The fact that its value, ( f4.7| ), is 
relatively large indicates that non-inclusive corrections to the shape function ( }4.5|) are 
important indeed. 



o Aleph 
□ Delphi 
o L3 
a SLD 

- NLL+power corrections 

- NLL 





(b) 



Figure 1: Heavy jet mass (a) and C— parameter (b) distributions at Q = M z with and 
without power corrections included. 

Having determined the parameters of the shape function, Eq. ( |4.7| ), at the reference 
energy scale Q = Mz, we can now apply the factorized expressions for the differential 
distributions, (fill), (|4.2|) and ( |4.3| ) with the same ansatz for the shape function ( |4.6|) to 
obtain the QCD predictions at different energy and compare them with the data. The 
combined plot for the p— and C— parameter distributions over the center-of-mass energy 
interval 35 GeV < Q < 189 GeV is shown in Fig. ||| a and b, respectively. Similar plot 
for the thrust distribution can be found in [nj. We observe that the theoretical curves 



reproduce the data over the whole interval of the shape variables including the end-point 
region. 



5. Moments of the event shapes 

Recently, the experimental data for the first few moments of various event shape distri- 
butions became available |]J . Their analysis indicates a presence of large hadronization 
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Figure 2: Comparison of the QCD predictions for the heavy jet mass (a) and 
C— parameter (b) distributions with the data at different center-of-mass energies (from 
bottom to top): Q/GeV = 35 , 44 , 91 , 133 , 161 , 172 , 183 , 189, based on the shape func- 
tion. 

corrections whose form deviates from IR renormalon models describing the nonpertur- 
bative corrections to the distributions as the shift of perturbative spectrum. 

Let us apply the obtained expressions for differential distributions to calculate the 
first two moments of the t— , C— and p— distributions defined as 

(e») = / dee 11 --, (n = 1 , 2) . (5.1) 
Jo °tot de 

Here, integration goes only over the part of the available phase space, < e < e max , 
corresponding to the three-particle final states, and it does not take into account the 
contribution of multi-jet final states, e > e max . Quantitative description of hadronization 
corrections to such final states is not available yet. Putting an upper limit on the value 
of the shape variable in fl5.1|) allows us to avoid the latter contribution and to replace 
the differential distribution da/de in fl5.1|) by the obtained expressions (|4.1| ), (|4.2|) and 
which are valid for < e < e max . 
Using general expression ( |2.10| ) one calculates the mean value of the event shape as 



(e) = (e) PT + ^ 



1 _ da PT (e n 

J- ^ max ' 



de 



+ 0(4) (5.2) 



Q 



where (...) PT = Jo™ de (...) da PT /de denotes averaging with respect to perturbative 
distribution and the scale (e) is defined as the first moment of the shape function, 
(s) = J deef(e). It is important to remember that the factorized expressions for the 
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differential distributions (|4.1|), ( |4.2| ) and fl4.3| ) are valid up to 0(1/ (Q 2 e)) — corrections 
which may modify the mean value (e) by 0(1/ Q 2 )— terms. The additional factor in 
front of (e) / Q takes into account that close to the edge of the three-particle phase space, 
e — ► e max , the perturbative distribution ( |2.3|) could take nonzero values 



d<? PT (e max ) a s (Q) 



de 



2n 



+ 0(a 



(5.3) 



This can be checked using the explicit expressions for perturbative distributions 
We find that da PT /de vanishes to one-loop order as e — > e max for the t— and p— variables 
while for the C— parameter it approaches a finite value. Finally, calculating the mean 
values (e) with respect to nonperturbative distributions f t (e) and f p (e,eQ) defined in 



Sp and (|4.4j ), respectively, we obtain 

(t) = (t) p . 

(p) = (P) 



pt + ^ + 0(1/Q 2 ) . 



Similarly, for the mean value of the C— parameter we get 



a s (Q) 
2vr 



5.73 + 0(a 



+ 0(1/Q 2 



(5.4) 
(5.5) 

(5.6) 



Here, large perturbative coefficient originates from (2.4) and it reduces a magnitude 
of the nonperturbative scale Ai by 11% at Q 



Mz- Relations (|5.4|) and (|5.5f) coincide 
with IR renormalon model predictions ^ ^j, O], while (|5.6| ) differs by perturbative 
a s (Q)— dependent "boundary" term. Nonperturbative Q— independent scale Ai is given 
by 

Ai = / de R / de L (e R + e L )f(e R ,e L ) = / deef t (e). (5.7) 



Substituting expression for the shape function, Eq. ( |4.6| ), one finds Ai = A(p(a,b) with 
ip given by 2F1— hypergeometric series. Using the values of the parameters (f4.7|) we find 



Ai = 1.22 GeV. 



(5.8) 



We would like to recall that this value depends on the factorization scale fi, Eq. (|3.12f ), 
and its \i— dependence is described by QCD evolution equation [l(|. Obviously, the 
value of Ai, and as a consequence 1/Q— corrections to the mean values (|5.4p, ( |5.5| ) and 
( p.6| ) are less sensitive to the choice of the parameters a, b and A as compared with 
nonperturbative corrections to the corresponding differential distributions. 

The comparison of the QCD predictions, (|5.4| ), ( |5.5|) and ( |5.6|) , with the data over 
the energy interval 35 GeV < Q < 189 GeV is shown in Fig. |3[ One should notice 
that (|5.4j ), (|5.5|) and ( ^.6|) describe the contribution of the two-jet configurations while 
experimental data take into account all possible final states. A good agreement observed 
in Fig. |] indicates that the contribution to the mean values of the final states with 
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Figure 3: Comparison of the QCD predictions to the mean values (t), (p) and (C) with 
the data. Dotted lines denote 0(a 2 ) — perturbative contribution, solid lines take into 
account power corrections given by Eqs. (|5.4|), (|5.5|) and (|5.6|). 



three and more jets as well as 0(l/(Q 2 e)) subleading corrections to the distributions 
are subdominant. Indeed, the dominant contribution to the moments ( |5.1| ) comes from 
the vicinity of peak of the differential distribution e n da/de. Using existing experimental 
data one can show [HJ that for n = 1 and n = 2 the position of the peak is located 
in the two-jet kinematical region while for higher n it moves towards larger e for which 
the integral (|5.1|) becomes very sensitive to the choice of the upper integration limit 
e max . We shall use this observation calculating the second moment of the event shape 
distributions. 

Let us apply ( |5.1| ) to calculate (e 2 ). Neglecting l/(Q 2 e) — corrections to the distribu- 
tion ( P-lOP we find after some algebra the following general expression 



(e 2 ) = (e 2 ) PT + 



+ 



(f!> 

Q 2 



2\ r 



Q 

1 Grr 



2(e) PT - e 



2 ^ PT (fimax ; 



de 



+ 



d(J pT (6 max ) 1 2 d 0" pT (6 max ) 



de 



+ -ei 



de 2 



(5.9) 



Similar to the mean value, (|5.2f) , the boundary terms vanish for the t— and p— variables 
while for the C— parameter they provide a sizeable contribution. Using the explicit 
expression for the shape functions, ([3.8|) and Q4.4Q, we calculate the scales (e 2 ) and take 
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Figure 4: Comparison of the QCD predictions for the second moments (t 2 ), (p 2 ) and 
(C 2 ) with the data. Dotted lines denote 0(a 2 s ) — perturbative contribution, solid lines 
take into account power corrections given by Eq. (|5.10D . 



into account the boundary terms Q2.4Q to obtaii 



Q T ' Q 2 



(P 2 > = (P 2 )P T + ^WP T + A2+ 4 Q2 2W) 

<C 2 > = (C\ T + || [2<0>„ - 



4.30 



+ 



9tt 2 A 2 

"q 2 



2tt 



11.46 



Here, the scale Ai was defined in ( p.7| ) and new scales A 2 and <5A 2 are given by 



A 2 = ({sr + e L f) , 6X 2 (Q) = ((e R - e L f) I I + 4 



pmax 



dp'p' 



da 



dp' 



(5.10) 



(5.11) 



where average is taken with respect to the shape function f(£R, El). The Q— dependence 
of the scale 5X 2 is attributed to perturbative prefactor depending on the single jet distri- 
bution, daY /dp, defined in flOP . Its origin was explained in Sect. 4. We find that the 



4 We are grateful to O. Biebel and S. Kluth for providing us 0(a 2 s )— expressions for the moments of 
the event shapes. 
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value of this factor varies from 2.19 at Q = 10 Gev to 1.85 at Q = 100 Gev. It is impor- 



tant to notice that ((er — El) 2 ) vanishes if one does not take into account non-inclusive 
corrections to the shape function (|4.5|) , 5/ n0 n-mci = 0. Using (|4.6|) and (|4.7|) one gets 

A 2 = 1.70GeV 2 , ((e r - E L f) = 0.14 GeV 2 . (5.12) 

It follows from ( |5.10|) that the boundary terms generate a sizable perturbative corrections 
to the second moment of the C— parameter distribution and diminish the magnitude of 
scales parameterizing l/Q— power corrections. Moreover, one finds from (|5.4|) , (|5.5|) , (|5.6| ) 
and ( J5.10 ) that the variance of the distribution, (e 2 ) — (e) 2 , does not receive l/Q— power 
corrections for the t— and p— variables while for the C— parameter the boundary terms 
produce a negative l/Q— correction 

(C 2 > - <C> 2 = (C\ T - (C)l T - 3.23 ^a s (Q) + 0(1/Q 2 ) . (5.13) 

The comparison of the QCD predictions (|5.10|) with the experimental data is shown in 
Fig. [|. We would like to recall that the obtained expressions for the moments do not 
take into account the contribution of multi-jet final states configurations and assume a 
smallness of l/(Q 2 e) — corrections to the distributions Q2.6Q . It is interesting to note that 
the last assumption is supported by the recent analysis of the power corrections to the 
first two moments of the thrust distribution in the single dressed gluon approximation 
| 20fl . This analysis is complimentary to our studies since it does not resum leading power 
corrections in the two-jet region and takes into account the contribution coming from 
the region e -> e max . 



6. Conclusions 

In this paper we have studied the power corrections to the thrust, t, heavy jet mass, p, 
and C— parameter distributions in the two-jet kinematical region. Our analysis was based 
on the observation |TtJ that perturbative and nonperturbative effects can be separated in 



the differential event shape distributions into calculable Sudakov resummed distribution, 
da-pT/de, and nonperturbative shape function, f(£R, £l), respectively. Each of them 
depends separately on the factorization scale fi but this dependence cancels in their 
product. The shape function describes the energy flow into two hemispheres in the final 
state. It does not depend on the center-of-mass energy Q as well as on the choice of the 
event shape variable e = t , p and C. 

We demonstrated that away from the end-point region, e 3> Aq C d/Q, nonperturba- 
tive corrections to the distributions have a simple form ( |2.8| ) with the leading l/Q— power 
correction parameterized by a single scale given by the first moment of the shape func- 
tion. In this region, to which all performed experimental analysis have been restricted 
so far, our predictions for the thrust and heavy mass distributions and their mean val- 
ues coincide with those of IR renormalon based models while for the C— parameters we 
find an additional sizeable perturbative contribution modifying the magnitude of the 
l/Q— power correction ( p.6| ). 
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In the end-point region, e ~ Aqcd/Q, the obtained factorized expressions for the 
distributions take into account power corrections of the form l/(Qe) n for arbitrary n. 
They are controlled by the shape function through (|2.11 ) and are sensitive to the choice 



of this function. Comparing the QCD predictions with the data we have chosen the 
simplest ansatz for the shape function ( }4.6|) which is consistent with general properties 
of nonperturbative distributions and includes nonzero correlations between energy flows 
into different hemispheres. Examining the dependence of the distributions on the corre- 
sponding parameter of the shape function we have observed that these correlations play 
an important role and are not negligible. 
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